Low-bandwidth control scheme for an oscillator stabilized Josephson qubit 
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We introduce a new Josephson junction circuit for which quantum operations are realized by 
low-bandwidth, nearly adiabatic magnetic-flux pulses. Coupling to the fundamental mode of a 
superconducting transmission line permits a stabilization of the rotation angle of the quantum 
operation against flux noise. A complete scheme for one-qubit rotations, and high- visibility Ramsey- 
fringe oscillations, is given. We show that high visibility depends on passing through a portal in 
the space of applied fluxes, where the width of the portal is proportional to the ramp-up rate of the 
flux pulse. 



Among the many candidates for the physical imple- 
mentation of a quantum computer, Josephson junction 
circuits have always been among the most promising. 
The quantum behavior of these circuits is readily tai- 
lorable by the choice of electrical topology and circuit 
parameters: there are various regions of this parameter 
space in which a coherent, controllable two-level quantum 
system, suitable for the realization of a qubit, is possi- 
ble. This same tailorability must also be exploited to 
avoid strong coupling to the environment and other de- 
cohering effects. The complexity of this optimization^] is 
such that many distinct Josephson circuit qubits are un- 
der active study, with successful single qubit control and 
two-qubit coupling achieved in a number of cases 0, Q 
But further improvements in qubit performance are un- 
questionably needed and continue to be sought. Recently, 
another important degree of freedom has been added to 
this search: strong, coherent coupling between a Joseph- 
son junction qubit and a quantized harmonic oscillator, 
realized by a, superconducting transmission line, has been 
achieved 0,0. This exciting discovery raises the question 
of how best such a coupled system is to be exploited for 
quantum information processing. 

In this Letter we report a new class of Josephson flux 
qubits that enter novel regions of the design space to 
achieve superior qubit performance. Here are its features: 
First, our qubit can be placed in a "frozen" state in which 
the barrier is very high, which makes resetting and mea- 
suring the qubit very reliable. Second, and more impor- 
tant, all qubit operations are realized by low-bandwidth, 
nearly adiabatic operations. Since the amount of envi- 
ronmental noise seen by the qubit is proportional to this 
bandwidth, we gain significantly by requiring only an ap- 
proximately 1 GHz control bandwidth rather than the 
many GHz that are necessary in other control schemes, 
which require the transmission of microwave radiation 
to the qubit. Finally, within our adiabatic operation 
scheme, the presence of a coupled quantum harmonic os- 
cillator plays a specific role in the qubit operation: the 
qubit rotation angle, which is equal to the dynamical 
phase difference accumulated between ground and ex- 
cited state energy surfaces, is stabilized by the adiabatic 
conversion of the qubit states to the ground and excited 



states of the quantum oscillator. 

The resulting qubit with its associated control trans- 
mission lines is complex, but the advantages gained by 
our operation strategy leads to a qubit for which a scale 
up to larger systems will eventually be more feasible. Our 
experimental realization of this qubit system is shown in 
Fig. 1 . The device consists of three Al / AI2 O3 / Al Joseph- 
son junctions grown using an in-house shadow mask pro- 
cess, arranged in a gradiometer pattern. This design is a 
modification of a qubit previously reportedQ. The body 
of the gradiometer is aluminum and consists of three 
loops. The lower loop is threaded by external flux $, 
the small left loop is threaded by the flux <E>c; the upper 
loop, the "pickup" loop, is inductively coupled to a high- 
Q niobium superconducting microstrip "pickup" trans- 
mission line. There are no wires attached to the qubit. 
Readout is done using a dc SQUID inductively coupled 
to the pickup transmission line. The qubit is operated at 
30 mK; however the effective electrical temperature from 
the circuits that drive the fluxes $ and $c* is about 1.3 
K. These circuits are not shown; however, these are cou- 
pled to the qubit with two additional superconducting 
microstrip transmission lines. The lines are shorted at 
the end that is not connected to the room temperature 
electronics. The length of the lines between the qubit 
and the shorted ends are such that at the desired fre- 
quency of operation of the qubit, 1.54 GHz, there will be 
a current node at the location of the qubit. Hence, when 
we operate the qubit at the degenerate point, which we 
can tune to have an operating frequency of 1.54 GHz, the 
contributions from both the low-frequency noise and the 
noise at the operating frequency arc minimized. 

Using our network graph theory @, we have obtained a 
comprehensive quantum mechanical model for this qubit. 
The system Hamiltonian is equivalent to that of a mas- 
sive particle in a four-dimensional potential. The particle 
mass is set by the system capacitances; the four degrees 
of freedom are the three Josephson junction phases plus 
the phase (the time integral of a voltage) across the ca- 
pacitance of the transmission line LC resonance. Other 
environmental degrees of freedom, e.g., those associated 
with the control transmission lines, are modelled as in 
6] using an oscillator bath. We calculate that the in- 
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FIG. 1: (a): Schematic layout of the qubit and the harmonic 
oscillator equivalent of the transmission line, (b): Picture 
of the qubit. The qubit is 500 /im long and each junction 
is 250 by 250 nm 2 in size. On the upper loop the pickup 
transmission line is shown, (c): Mask art of the three shadow 
type Josephson junctions (blow up of dashed box in (b)). 



trinsic quantum lifetimes associated with this bath, for 
the choices of control parameters described below, and 
for T = 100 mK, are T 2 = 600 nsec and T x = 300 nsec, 
which are more than long enough for our adiabatic con- 
trol scheme to successfully execute quantum operations. 
In practice the coherence of our qubit is further decreased 
by extrinsic noise in the control electronics, an example 
of which will be analyzed below. 

While the resulting potential function is complex, there 
is a region of two-dimensional parameter space of applied 
fluxes $ and <&c (see Fig. 1) in which the potential has a 
simple, double- well structure. The lowest energy states of 
this double well form the qubit; the qubit is also linearly 
coupled to the transmission line harmonic oscillator, as 
we will describe shortly. 

In the vicinity of these double-well minima, the three 
degrees of freedom of the qubit can be linearly trans- 
formed into two "fast" degrees of freedom transverse 
to the double-well axis and one slow longitudinal phase 
which we will call 5. Using a Born-Oppcnheimer 
approximation^ for the fast degrees of freedom, we 
have numerically obtained an effective one-dimensional 
double- well potential Vq(6) for the qubit. While we use 
the full numerical form in our calculations, it is useful for 
describing features of this potential to represent it using 
the following simple analytic anharmonic formj^]: 



Vq{6, t) = -h 2 (t)6 2 + h 4 6 4 + a(t)S. 



(1) 



$c in time. In our devices, to good approximation, $ 
controls the double-well asymmetry coefficient a and <&c 
varies the double- well barrier height h% . We find that h% 
can be either positive or negative; that is, the double- well 
structure can be transformed into a single well. 

Eq. (JIJ is only one part of the description of the sys- 
tem; coupling to the harmonic-oscillator phase (p must 
also be included. To good approximation, the full Hamil- 
tonian of the system can be written 
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We indicate here an explicit dependence on time t, as 
our qubit control technique involves the pulsing of $ and 



Here the first two terms are the qubit Hamiltonian, the 
next two terms are the harmonic oscillator Hamiltonian, 
and the last term is the linear coupling between the 
two. Q and q, the qubit and oscillator charge oper- 
ators, act as momentum operators in the Schrodinger 
equation: [ip, q] = [<5, Q] = i. Cq is the capacitance as- 
sociated with the Josephson junctions; in our circuit all 
three have Cq « 50fF (there is substantial excess ca- 
pacitance in parallel with the intrinsic capacitance of the 
oxide junctions). Ct ~ IfF and Lt ~ lOnH are the 
effective transmission line capacitance and inductance, 
Lq sa QAQpH is an effective qubit inductance (actually a 
complicated function of the full inductance matrix of the 
qubit circuit 6]), M sa 50pH is the mutual inductance be- 
tween the qubit and the transmission line, and <E>o = h/2e 
is the superconducting flux quantum. In other studies^] 
the physics of this situation has been described using a 
Jaynes-Cummings Hamiltonian; we find that retaining 
the first-quantized form of Eq. J5J has proved helpful in 
analyzing the details of our control scheme. 

In our device, the barrier height I12 in the model po- 
tential Eq. 0}, at low values of the control flux 
can easily be made so high that the quantum tunnelling 
rate A can be made very small (S/A being many min- 
utes at least). In this regime the qubit state is essentially 
frozen, with two degenerate eigenstates \L) and \R) for 
"left" and "right" well, see Fig. 3, point A. The device is 
hysteretic in this regime, so it can be stably set in either 
\L) or \R) by a suitable quasistatic sweeping of flux (f>. In 
addition, the state is easily read out by a sensing SQUID 
inductively coupled to the transmission line of Fig. 1. 

So, for $c = $ = 0, corresponding to a, A w 0, the 
qubit is in a "resting condition" in which its dynamics are 
frozen. Starting from this resting condition, conceptually 
there is a simple prescription for adiabatically applying 
single-qubit quantum gates to the qubit. A U Z rotation" 
(of the Bloch sphere defined by the \L, R) basis) is simply 
obtained by adiabatically varying a in Eq. (J) (by pulsing 
the main-loop flux <3>) with A held equal to zero. In the 
adiabatic limit, the Z rotation angle is given by <f> = 
2 / dte(t)/h, where 2e = ((L\Vq\L) - (R\V Q \R)) is the 
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difference in depths of the two wells of Vq . 

The ll X rotation" scheme will require more discussion, 
as it involves a stabilizing effect due to the harmonic os- 
cillator degree of freedom in Eq. J2J • I£ starting from the 
resting state, the control flux <j>c is pulsed, then the e de- 
fined in the last paragraph will remain zero and hi will be 
varied. Since ideally the double-well potential Eq. JTJ is 
symmetric throughout this pulse, the lowest-lying eigen- 
states of this potential will always be symmetric (\S)) 
and antisymmetric (\A)) with respect to the qubit coor- 
dinate S. In the resting state \S, A) = (\L) ± \R))/V2; 
these basis states lie at right angles to the \L, R) states on 
the Bloch sphere, so that adiabatic evolution along this 
coordinate and back results in an X rotation by angle 
9 = J dtA(t)/h, where A(t) = £i($ c (i)) - £ ($c(t)), 
and are the first two instantaneous eigenenergies (see 
Fig. 2) of the qubit Hamiltonian Eq. (J2J. 
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FIG. 2: Instantaneous eigenvalues of H, Eq. as a function 
of control flux <3?c, for main flux $ ~ 0. The exponentially 
rising tunnel splitting A, and the $c-independent oscillator 
frequency u>t, are indicated. At A, the "resting state", the 
double- well barrier is very high and the two lowest eigenstates 
are nearly degenerate; here the resting state is not quite sym- 
metric, with e w 15MHz. For this nonideal case successful 
quantum operation requires that the system pass through the 
"portal" as described in the text. At point B, beyond the 
portal, quantum tunneling dominates the dynamics and the 
system is effectively symmetric. Beyond B the qubit comes 
into resonance with the transmission line oscillator, with a 
splitting of 2g ~ 220MHz. At point C, the lowest lying states 
have purely harmonic oscillator character. 

In Fig. 3 we see that the first two eigenfunctions in 
the potential of Eq. J2J change form in going from point 
A to point B, from \L,R) to \S, A). This rapid change 
is actually desirable for the operation of the qubit; nu- 




FIG. 3: Color density plot of first three first energy eigenstates 
of H at points A, B, and C of Fig. 2. 5 is the horizontal axis, 
and (p the vertical axis (see Eq. For A the first two 

eigenstates, localized in the left and right wells, are not quite 
degenerate because e / (see Fig. 2). In going from A 
to B, on account of strong quantum tunneling, the first two 
eigenstates change rapidly in character, becoming symmetric 
and antisymmetric combinations of the well states. For both 
A and B, the third state has one quantum of excitation of the 
harmonic oscillator. The states at C have purely harmonic- 
oscillator character. In going from A to C, the steady decrease 
of the well-well separation is evident. 



merical integrations of the time-dependent Schrodinger 
equation show that we obtain high-fidelity quantum op- 
erations if $c(t) is pulsed with an overall rise time of 15 
nsec, with appropriate pulse shaping (see the discussion 
of visibility-limiting mechanisms below). 

Without the coupled harmonic oscillator, this X- 
rotation scheme suffers from the problem that A, and 
thus 9, is exponentially sensitive to the value of <&c- This 
exponential sensitivity to the barrier height is evident in 
Fig. 2 for values of $c near B. However, the presence of 
the oscillator degree of freedom changes several features 
of the spectrum in a dramatic way. The essentially hor- 
izontal lines in Fig. 2 are the equally spaced oscillator 
energy levels, which are insensitive to the flux applied 
to the qubit. But because of the mutual inductance M 
in the last term of Eq. J2J), there is an avoided crossing 
between the oscillator and qubit eigenlevels. For the pa- 
rameters of the device in Fig. 1, the "vacuum Rabi split- 
ting" 2g (in the notation of the Jaynes-Cummings model, 
see 0|) indicated in Fig. 2 is about 220 MHz. For the 15 
nsec ramp-up time mentioned above, the time-evolution 
remains accurately adiabatic through this anticrossing as 
well. Thus, after ramping up &c to point C in Fig. 2, 
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the system remains in a superposition of just two states, 
the ground state and the first excited state which now 
has almost entirely the character of a single quantum of 
excitation in the harmonic oscillator. A in this regime 
is almost completely independent of $C) meaning that 
the execution of the X rotation in this regime is almost 
entirely insensitive to noise in $c- 

While not of direct relevance to small-angle X rota- 
tions that will be desired for quantum gates, this insen- 
sitivity is very valuable for performing the analog of a 
Ramsey- fringe experiment to characterize this qubit 
The protocol for this experiment is as follows: First, pre- 
pare the system in the state \L), which in the ideal case 
(e = 0) is an equal superposition of the energy eigen- 
states \S) and \A). Second, ramp up the control flux 
$c adiabatically, then hold it constant at a value <&hoid 
past the anticrossing with the harmonic oscillator state 
(i.e., to point C in Fig. 2). After hold time T, ramp 
down the control flux to the resting state, then measure 
the qubit in the \L,R) basis. The measured L probabil- 
ity should show Larmor oscillations of the form coslutT, 
where lot = 

We believe that the inaccuracy in the setting of the 
resting state, i.e., the inability to set e exactly equal to 
zero, is likely to be the most important mechanism for 
reducing visibility in our system. If the resting state is 
not exactly symmetric, then, since A is exponentially 
small at point A, it is likely that A << e. In this case 
\L) is not an equal superposition of energy eigenstates, 
but is an eigenstate itself, as Fig. 2a shows. Thus, a 
completely adiabatic evolution will keep the system at 
all times in an instantaneous eigenstate, and no Larmor 
oscillations will take place - the visibility will be zero. 

The resolution to this dilemma is that, for the realistic 
case in which f / 0, the time evolution in the Ram- 
sey experiment we have discussed is at a crucial moment 
actually nonadiabatic, in a way that permits the experi- 
ment to succeed as described. There is a "portal" in the 
e-A parameter space, which, if passed through at a finite 
rate, results in a successful X rotation (Ramsey-fringe 
visibility near 100%). The idea is that since A rises ex- 
ponentially, A(i) = Aoe 4 / T , there is a brief interval of 
time during which A(i) « e; this sets the location of the 
portal in the parameter space. This is the critical time 
because here the form of the eigenstates of H changes 
quite rapidly from nearly \L, R) to very nearly \S, A), as 
Fig. 3 shows. The evolution here will be nonadiabatic 
if er/Ti is not too large. That is, —1 < er/Ti < 1 ap- 
proximately defines the width of the portal. If the time 
evolution passes through the portal, then nearly half of 
the state amplitude is transferred from the ground state 
to the excited state, resulting in essentially 100% visibil- 
ity. 

This picture has been confirmed in detail by numer- 
ical simulations, and also emerges from an exact solu- 
tion of the two-state problem in which e is time inde- 



pendent and A increases exponentially [9j; then the ex- 
act solution of the time-dependent Schrodinger equation 
is proportional to c+(t)\L) +ic-(t)\R), where c±(t) = 
exp(t/2r)J ± i/2 + i eT /R[A r/?iexp(t/r)], and J v (x) is a 
Bessel function. With this exact solution the visibility 
for finite e can be calculated to be sech 2 (7rer/2?i). Our 
numerical simulations confirm this dependence: given the 
existing precision of pulsed magnetic-flux control (a few 
/i<i>o)j a nonzero e in the MHz range will typically be 
present. This implies that the portal will typically be 
located near <j>c = 0.35<i>o as in Fig. 2. We numerically 
find a Ramsey- fringe visibility in excess of 90% if <I>c if) 
is initially ramped rapidly through the portal (from A to 
B in 0.8 nsec), and then brought slowly up to its final 
value (from B to C in 15 nsec). The necessary time de- 
pendent flux control is demanding, but can be achieved 
in our lab. 

We believe that our qubit and its control schemes show 
good promise for scalability to larger systems. The rel- 
atively large size of the main loops (Fig. 1) means that 
establishing strong inductive couplings between qubits 
is straightforward, permitting fast two-qubit gates to be 
done. As we have already demonstrated, its large size 
also makes reliable, and potentially fast, readout of the 
qubit state quite simple. The strong coupling to a trans- 
mission line degree of freedom permits X rotations to be 
done that are very insensitive to flux noise, and, as other 
workers have recently shown, the presence of a control- 
lable harmonic-oscillator quantum also offers new possi- 
bilities for quantum computing architectures, including 
easy motion of qubits and long-distance coupling of sta- 
tionary qubits. Our control schemes now considerably 
enrich this picture. 
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